Role of seeding the cavity of a two-photon correlated emission laser with thermal light 
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A study of the evolution of two-mode squeezing, entanglement and intensity of the cavity radiation 
of a two-photon correlated emission laser initially seeded with a thermal light is presented. The 
dependence of the degree of two-mode squeezing and entanglement on the intensity of the thermal 
light and time is found to have more or less a similar nature, although the actual values differ 
specially in the early stages of the process and when the atoms are initially prepared in nearly 50:50 
probability to be in the upper and lower energy levels. Particularly, seeding the cavity turns out 
to spoil the nonclassical features significantly in the vicinity of t = 0. It is also shown that the 
mean photon number in a wider time span has a dip when mode b is seeded, but a peak when 
mode o is seeded. Moreover, this study asserts that the effect of the seeded light on the nonclassical 
features and intensity of the cavity radiation is eroded with time by the pertinent emission-absorption 
mechanism which can be taken as an encouraging sign in practical utilization of this quantum system 
as a source of bright entangled light. 



PACS numbers: 42.50.-p, 42.50.Ar, 42.50.Gy, 42.50.Lc 



I. INTRODUCTION 

A three-level cascade laser has received a considerable 
attention over the years in connection with the strong 
correlation between the modes of the generated radiation 
that leads to a substantial degree of nonclassical features 
[H4llj]. One of the various possible schemes, in this re- 
gard, is a correlated emission laser in which the atomic 
coherence is induced by initially preparing the atoms in 
the coherent superposition of the atomic energy levels 
between which a direct spontaneous transition is elec- 
tric dipole forbidden (description of the model can be 
found for instance in [H, 0, EH)- This quantum system 
has been thoroughly studied elsewhere for different sit- 
uations when the cavity is assumed to be initially in a 
vacuum state [J, S, E3, III ■ 

It is observed in these works that a strong robust en- 
tangled light can be produced by manipulating the rate at 
which the atoms are injected into the cavity and the ini- 
tially prepared coherence. In addition, the mean photon 
number was found to be higher where the degree of en- 
tanglement is larger (l3| . However, the assumption per- 
taining to no radiation initially in the cavity appears to 
pose a serious challenge when weighed solely from tech- 
nical point of view since it is practically impossible to 
build a cavity entirely in the vacuum state. Even if one 
is determined to do so, it requires a significant effort and 
resource to cool it at that level. Undoubtedly, this recurs 
additional cost leaving the emerging technical difficulty 
aside. 

To circumvent the prevailing predicament somehow, 
recently there are efforts in assuming the cavity to con- 
tain a coherent light with certain number of photons 
when the atoms are taken to be pumped externally with 
another coherent light [1, [l4| • It was predicted that the 
coherent light increases the mean photon number and en- 
hances the dependence of the entanglement on the phase 



fluctuations. At this juncture, there is one outstanding 
issue to be addressed in relation to having three coher- 
ent lights (coherent light with different frequency) in the 
cavity at the same time. 

From naive understanding of atom-radiation interac- 
tion, it would be reasonable to expect the atoms to re- 
spond to the coherent lights in the same manner. Conse- 
quently, as long as there is a sufficient amount of resonant 
coherent light in the cavity, there is a good chance of initi- 
ating a two-step pumping process. This provides another 
root for inducing a coherent superposition between the 
lower and the upper energy levels that amounts to estab- 
lishing a similar effect as the external coherent radiation. 
In light of this, the general tendency of designating the 
observed result to the pumping radiation and the initial 
coherent radiation in the cavity seems to beg for a more 
serious indepth consideration than the mere inclusion of 
the coherent light as initial state in the characteristic 
function can warranty. 

One of the most realistic proposals to study the effect 
of the light in the cavity at the beginning on the gener- 
ated radiation may be heating the cavity to certain par- 
ticular temprature using ordinary light source and then 
inject the initially prepared atoms to it. This can be done 
by placing the cavity in a bath of ordinary light allowed 
to leak into it via the coupler mirrors and wait until ther- 
mal equilibrium is reached or put the light source in the 
cavity until a wall temprature of interest has reached and 
then turn it off prior to the injection of the atoms. The 
former scenario is the usual consideration of the thermal 
reservoir treated elsewhere .12]. In this case, the ther- 
mal light is found to increase the mean photon number, 
but degrade the two- mode squeezing and entanglement. 
However, the analysis of the latter is yet to be addressed. 

With this understanding, in this communication, the 
effect of the thermal light initially seeded in the cav- 
ity as in the second option on the statistical properties 
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and quantum features of the cavity radiation is studied 
where the thermal light is conveniently represented by 
the chaotic state. Although the coupling of the quan- 
tum system to the environment invariably disparages the 
entanglement and other quantum features, the effect of 
decoherence due to thermal light is not taken into con- 
sideration. Nevertheless, the environment is assumed to 
be made of vacuum modes whose fluctuations would be 
included in the usual manner. 

Since the approach followed is presumed to be adapt- 
able for arbitrary reservoir with a minimum effort, the 
thermal light outside the cavity has not been given at- 
tention to deter the confusion. Basically, the cavity con- 
structed in an open space at room temprature and then 
placed in an experimental setup can be considered as the 
cavity that contains chaotic light. In addition to this, 
when it becomes practically impossible to pump out all 
radiation from the cavity, so that it would be in a vac- 
uum state, the remnant light can also be treated as a 
chaotic light. Whatever the case may be, studying the 
effect of light in the cavity before the injection of atoms 
is envisioned to be of a paramount importance in utiliz- 
ing the generated radiation for particular purposes like 
teleportation and quantum information processing as an- 
ticipated. 

To achieve the intended goal, the procedure by which 
the chaotic light in a cavity can be integrated with the 
three-level atomic system would be formulated. To do so, 
the general characteristics of thermal light is outlined in 
the way that it can be salvaged for present work. Then 
the equations of evolution for the three-level laser for an 
arbitrary initial radiation would be obtained where the 
contribution of the thermal light enters via this initial 
condition. In the process, it is assumed that the atoms 
are unable to distinguish among cavity radiations other 
than their frequency whereby the thermal light can also 
be absorbed by the atoms. In light of this, the effects 
of initially heating the cavity (thermal seeding) on the 
achievable two-mode squeezing, entanglement and mean 
photon number is investigated by making use of the cor- 
relation properties of the chaotic light. 



II. SKETCHY DESCRIPTIONS OF THE 
INVOLVED SYSTEMS 



A. Chaotic light 

It is a well established fact that the radiation gener- 
ated by conventional sources of light can be taken as the 
best example of chaotic light. For instance, the light 
produced by a laser machine operating below threshold 
exhibits chaotic nature. For the sake of convenience, let 
us consider the thermal light in a cavity maintained at 
equilibrium with the walls at temperature T . The den- 
sity operator p for the thermal light can be related to the 



entropy by von Neumann relation [151 ] 

S = -k b Tr(p Inp), (1) 

where Kb is the Boltzmann constant, Tr stands for trace 
operation and p is the density operator. 

On the basis that at equilibrium entropy should be 
maximum and Tr(p) — 1, the density operator for the 
chaotic light is related to the temprature by 

-h/k b t 

P = : , (2) 

Tr{e- H /*B T ) 

where H = uia'a is the Hamiltonian of the radiation field 
when the zero point energy is shifted by w and h = 1 for 
convenience. 

The same density operator can be described in the 
number state applying the completeness relation for num- 
ber state / = X)^°=o l n )( n l an d the orthonormality con- 
dition (n\m) = S nm as 



P = 



£.71 



(1 + n) 1+n 



(3) 



where n is the mean photon number. Since the chaotic 
light does not exhibit correlation, for two-mode light, the 
density operator can be treated as separable; 
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(4) 



where h a and n b are the mean photon numbers corre- 
sponding to the respective modes. 

On the other hand, the number state can be describ 
able in terms of the vacuum state as 

, x (a + ) n 



|o>, 



(5) 



which can also be expressed in terms of the coher- 
ent state using the power series expansion of the form 

e — Z^n=0 



as 



oo 

E 

n=0 



(6) 



Applying the density operator (0| along with the rela- 
tion of number state with coherent state © and defini- 
tion of the mean value (O) = Tr(Op), for any operator 
O, it is possible to verify for two- mode chaotic light that 

(a*a) = n a , (7) 
(/3*/3) = n b , (8) 

(a(3) = (a 2 ) = (/3 2 ) = (a) = (p) = 0. (9) 

If it is found necessary, the mean photon number of 
each mode can be related to the frequency of the corre- 
sponding mode and temprature of the wall by 

' (10) 



i/K B T _ I 



with i = a,b. 



3 



B. Two-photon correlated emission laser 

A two-photon correlated emission laser can be inter- 
preted as the amplification of simultaneously emitted 
photons with different frequency that are correlated due 
to the coherence superposition induced by preparing the 
atoms in the upper and lower energy levels denoted by 
| a) and |c), respectively. In establishing the lasing mech- 
anism, the atoms are presumed to be injected with con- 
stant rate r a into the cavity and left after successfully 
decayed to other energy levels. 

In general, interaction of a nondegenerate three-level 
cascade atom with a two-mode cavity radiation can be 
described in the rotating-wave approximation and the 
interaction picture by the Hamiltonian of the form 

H AR = ig[a\a)(b\ - \b)(a\tf - \c)(b\tf + b\b)(c\], (11) 

where g is the coupling constant taken to be the same 
for both transitions, a (b) is the annihilation operator for 
one of the cavity radiations and \b) is the intermediate 
energy level. 

The atoms are initially prepared in arbitrary coher- 
ent superposition of energy levels between which a di- 
rect spontaneous transition is electric dipole forbidden, 
that is, the initial state of the atom is chosen to be 
IMO)) = C a (0)\a) + C c {0)\c), where C o (0) and C c (0) 
are probability amplitudes for the atoms to be initially 
in the upper and lower energy levels. The corresponding 
initial density operator would be 



p A (0) = pW\a)(a\ + p%>\a)(c\ + p^\c)(a\ + p™\c)(c\, 

(12) 



JO) 



(0)1 



where p^a 



and p^cc are the corresponding populations 
and p y ac represents the initial atomic coherence. 

With this condition, applying the linear and adiabatic 
approximation schemes in the good cavity limit, the time 
evolution of the reduced density operator that accounts 
for the contribution of the cavity radiation due to the 
atoms is found to be [la] 



„(o) 



III. EQUATIONS OF EVOLUTION OF THE 
COMBINED SYSTEM 



In order to study the dynamics of the cavity radia- 
tion of the combined system, it is necessary to obtain the 
associated equations of evolution. To this end, the con- 
tribution of the initial chaotic light in the cavity and the 
two-mode vacuum reservoir towards the master equation 
is sought for. 

It is a well established fact that the time evolution 
of the reduced density operator for the cavity radiation 
coupled to a reservoir has, in the Born approximation 
1171 . the form 



dp(t) 
dt 



= -i[H St m]-i[(HsR(t)) R , p(0)] 



[(H SR (t)) R , [H s (t'), p(t')]]dt' 



- / Tr R [H SR (t), [H SR (t'), p{t')R}]dt\ (14) 
Jo 

where S and R refer to the system and reservoir variables 
and p(0) represents the radiation initially in the cavity. 

Assuming that the already heated cavity is placed in a 
vacuum (environment modes with sufficiently small pho- 
ton number to disturb the analysis) and the injection of 
the atoms begins right after that, one can designate p(0) 
as a density operator for the chaotic light in the cavity. 
Furthermore, the interaction of a two-mode cavity radi- 
ation with a two-mode reservoir can be described in the 
interaction picture by the Hamiltonian 



Hsn(t) = i^A fe [a t a fe e 1 



+ &tfc fce i("o-^)t _ fat -i(u, -u k )tl 



(15) 



where ojq = u)a + u)b , with u a and Wb stand for the frequen- 
cies of the cavity radiations, (&&, bk) are the annihilation 
operators, Wk is the frequency and is the coupling 
constant for the k th mode of the reservoir. 
In view of Eq. (|15p , one can write 



(Hs R (t)) R = i Afc[a t (afc)fle 



i(uo— u>h)t 



dp A(l - rj) 



dt 4 
L - 

~4 



2a) pa — da) p — pad) 



+ + [2b0 -tfbp-ptfb] 



Ay/1 4 ^ [paW - 2atpSt + aWp 
2b pa + abp + pab\ , 



(13) 



a(al)« e " 4( "°^ fc)t + & t <&fe)«e i(w °- a " t)i 



-b{bl) R e- i ( ua - Uk ) t ]. 
For a two-mode vacuum reservoir, 



(afe)ij = {bk) R = 0, 



which leads to 



(16) 



(17) 



2r 2 

where A = ra .f is the linear gain coefficient, 77 — 



(0) 



Pec 



p^aa is the population inversion and 7 is the atomic decay 
rate. 



SR/R 



= 0. 



(18) 



This entails that the thermal light in the cavity does 
not directly contribute to the master equation. Hence 
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solving the remaining terms following the standard ap- 
proach [nj yields 



dp k 
~dt ~ 2 



\2dpa) — a) dp — pa) a] 



A(l-rj) rri A+AA AA+A AAA+n 
H [2a 1 pa — aa' p — paa' \ 

+ - [A(l + + 2k] [2bpb 1 - b r bp - 0b] 

+ ViEZ [^tgt _ 2a tp t + atfotp 
— 2o/5a + abp + pab] , 



(19) 



where k is the cavity damping constant. 

Employing the master equation (1191) and following a 
straightforward algebra, it is possible to obtain [l], 0] 

a(t) = A + (t)a(0) + B+(t)P*(0) + C+{t) + D+(t), (20) 



f3(t) = A_(t)/3(0) + B_(t)a*(0) + CL(t) + D_(i), (21) 
where 

A±(t) = ^ [fo T lK^ 4 + (ry ± lje"**] , (22) 



B±(«) = ± 



y/l — rj" r : L_ 



2,/ 



e 2 - e 2 



cw-ijfK,-!: 



e 2 



(i-t') 



+ (» ? + l)e-*( t - t ')] /„(«')< 



•j + Ai, 



(*-*') 



+ (ij-l)e-t( t - t ')]A(f)df, 



£>+(*) = 



n2 /■* 



e 2 



2r? 

■(*"*')] /*(f)df 



(23) 



(24) 



(25) 



(26) 



= (fa(tVa(t)} = 0. (30) 

It is worth noting that a(0) and /3(0) along with their 
corresponding complex conjugates represent the thermal 
light. It is clearly seen from Eqs. (|2U| and (|21l) that the 
thermal light enter into the evolution of the system as an 
initial condition as anticipated. But if the source of the 
chaotic light remains in the cavity throughout the lasing 
operation, the effect of the thermal light should also be 
included in the correlations of the noise forces. 



IV. EVOLUTION OF THE TWO-MODE 
SQUEEZING AND MEAN PHOTON NUMBER 

A. Two-mode squeezing 

Without lose of generality, a two-mode cavity radiation 
can be described by an operator 



1 



c= T2 {a + b). 



(31) 



Particularly, the squeezing properties of the cavity radi- 
ation can be studied applying the quadrature operators 
defined by 



c + = c) + c 



and 



c_ = i{& — c), 
with the commutation relation 

\c, , c 1 = 2i. 



(32) 
(33) 
(34) 



Applying Eqs. ((3TJ), ([32]), (J33J) and flMJ), the quadrature 
variances in terms of the c-number variables associated 
with the normal ordering are found to be 



Ac| = 1 + 1 [2(a*a) ± (a* 2 ) ± (a 2 ) 



(a*) 2 - (a) 2 



± (f3* 2 ) ± (f3 2 ) + 2(0*0) - (n 2 - W 2 ] ± (a*P*) 

+ (a*P) + (a/3*) ± (ap) T (a*) (a) - (a*)(/3*) 

T (a*)(0) T (a)(0*) - (a)(0) T (0*)(0). (35) 



D-{t) 



2 /■* 



e 2 



27? 



e 2 



1 (*-*') 



in which 



(27) 



(/-(*')/:(*)> = ^r-^*(*-* / ). (28) 



The various correlations in Eq. (|35l can be evaluated 
applying Eqs. (|20| and ([2l"j). Thus on the basis of Eqs. 
(|20| and ((2T|) , and assuming the cavity modes to be ini- 
tially in a two-mode thermal light, one can readily verify 
that 

(a(t)) = (/?(*)) = 0- (36) 
Moreover, upon taking Eq. (|20|) into account, 



(h(t')fa(t)) = Ay/ \^ 6(t-t>), (29) 



(a*(t)a(t)) =A+(t)(a*(t)a(0)) + B+(t)(a*(t)0*(O)) 

+ (a*(t)C + (t))+(a*(t)D + (t)). (37) 
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In line with this, using the complex conjugate of Eq. (l2Cfl) , 
it is possible to see that 

(a*(t)a(0)) = A+(t){a*(0)a(0)) + B + (t)(a(0)(3(0)) 

+ (a(Q)C* + (t)) + (a(0)D* + (t)). (38) 

Based on the fact that the noise force at time t does 
not affect the cavity mode variables at earlier times and 
taking the cavity modes to be initially in a chaotic state, 
one gets 

(a*(t)a{0)) =A+{t)n a . (39) 
It is also possible to verify in a similar manner that 

(a*{t)j3*(p))=B+(t)n b , (40) 

(a*(t)C + (t)) = (C* + (t)C + (t)) + (D* + (t)C + (t)), (41) 

(a*(t)D + (t)) = (C* + {t)D + (t)) + (D* + (t)D + (t)). (42) 

Hence upon substituting Eqs. ([39"]). (|40j) . (T4T1) and (|42|) 
into ([57)1 . one obtains 

(a*(t)a(t)) = (C* + (t)C + (t)) + (D* + (t)D + (t)) 
+ (C* + (t)D + (t)) + (D* + (t)C + (t)) 
+ A\{t)n a + B\{t)n b . (43) 

Next the remaining correlation functions in Eq. ([15)1 
would be evaluated. To this end, it is possible to write 
having Eq. flU} that 

(C* + (t)C + (t)) = -^f [\(r]- i )e -(»+^)(t-t')/2 
47 7 Jo Jo 

+ + l)e- K ( t -*')/ 2 ][(ry - i) e -(«+^)(*-*")/2 
+ (r, + l^-^/^Unf^dt'dt", (44) 
from which follows using Eq. 



(Cl(t)C + (t)) 



(r ! -l) 2 (l-e'^ +A ^ t ) 
16(k + Arf) 



+ 



(r, + l) 2 (l 



16k 



(77 2 -l)(l-e-( 2K + A ")*/ 2 ) 



(45) 



4(2k + Arf) 

It is not difficult to demonstrate in a similar manner that 
(D* + (t)D+(t)) = 0, (46) 

(C* + (t)D + (t)) = (D* + (t)C + (t)) 
_ A(l - rf) \( V - 1)(1 - e -(«+^)*) 



(tj + 1)(1 



32(k + A??) 



(1 



-(2k+Ajj)*/2 



32k 



8(2k + A//) 



(47) 



Therefore, in view of Eqs. J22]), ([23]), d43|), gS]), (|46j) 
and (|47[) . one can readily establishes 



<"*<»° <, »-^T|y (1 -^ , " +J "" ) 



(l-e' 



-(2 K +A)j)t/2^ 



2t](2k + Arf) ' 

"" f (r? - lfe"^^' + (tj + 1)V 
n b {l-rf) 



4?7 2 

+2(77 2 - i) e -( 2 «+^)*/2 



4?7 2 



g -(K+A»7)t + e -«i _ 2e -(2K+A J) )t/l 



(48) 



C9*(*)|8(t)) 



A(l - v 2 ) 



+ 



Ar)(n + Arf) 
A(l - V 2 ) 

2i](2k + Arf) 

fib 



n _ e -(«+^-»7)*) 

_ e -(2K+A»7)t/2-j 



(r ? + i) 2 e -(«+^) t + _ l)- e 



2„-Kt 



4?7 2 

+2(t7 2 - i) e -(a«+Au)t/2 



n a (l-7? 2 ) 
4^2 



x e -(K+^)« _|_ g-Kt _ 2e -(2K+Ar,)t/2 ^ Mg\ 



<a(t)j9(t)) 



A(l- V )y/1 -I) 2 _ e _( re+ A,) t) 

4?7(k + Arf) ' 



+ 



Ay/1- 


- ?7 2 


2t](2k + 


A77) 


n a \/l - 


r; 2 


Ami 2 




(r? + l)e" 


-Kt 




rf r 






(ry-l)e- 


-Kt 



(r] - i) e -("+A?)t 

2e -(2«+A^)t/2" 

(»? + i) e -(»+^)* 



(50) 



(a 2 W) = (/3 2 (i)) = (a*(t)/3W)=0. (51) 



On account of Eqs. (|36)) and (jBljl . the quadrature vari- 
ances reduce to 

Ac| = 1 + {a*(t)a(t)) + ((3*(t)(3(t)) ± 2(a(t)/3(t)}. 

(52) 

In order to investigate the time development of the two- 
mode squeezing, applying Eqs. (|ig|l. (gSJ), (I5U)) and ([ST]) . 
Ac 2 , is plotted against time for k = 0.5, r\ = 0.2 and 
various values of A given n a and ri&. The value of 77 is 
fixed based on earlier study at steady state in which the 
degree of squeezing is found to be the highest at 77 = 0.2 
for A = 10 [l| and recent report for arbitrary time 
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A=20 
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time 
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<n> =0.01 
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<n> =0.50 

b 

<n> =1.50 



FIG. 1: Plots of the minus quadrature variance (Ac 2 .) of the 
cavity radiation for n — 0.5, n a — hb — 0, r\ — 0.2 and 
different values of A. 



FIG. 3: Plots of the minus quadrature variance (Act) of the 
cavity radiation for k = 0.5, n a — 0, r\ — 0.2, A = 10 and 
different values of fib- 



2.5 


r 


<n>. 


=0.01 




I 


<n> 


=0.50 




I 
I 


<n> 


=1.50 


2 


i 







1.5 



0.5 



\ \ 



4 6 

time 




FIG. 2: Plots of the minus quadrature variance (Ac?.) of the 
cavity radiation for « = 0.5, fib = 0, r\ — 0.2, A = 10 and 
different values of n a ■ 



FIG. 4: Plots of the minus quadrature variance (Ac 2 .) of the 
cavity radiation for k = 0.5, r\ — 0.2, A = 10 and different 
values of n a = n&. 



It is clearly shown in Fig. Q] that the degree of two- 
mode squeezing increases with the rate at which the 
atoms are injected into the cavity and time. This en- 
tails that with time more atoms get chance to traverse 
the cavity whereby the nonclassical correlation acquires 
more strength. This type of reasoning may give sense if 
one compares the result obtained here with the earlier 
studies at steady state [if. It is also good to note that 
the degree of two-mode squeezing has peak value when 
i] = and there is an external pumping radiation [l8j . 
It has been argued elsewhere that such property of the 
two-mode squeezing is related to the thermal fluctuations 
resulting from vibrations of the atoms on the walls of the 
cavity. From reports at steady state, it is not difficult to 
observe that the mean photon number at rj = is sig- 



nificantly larger than at 77 = 0.2. In light of this, it is 
possible to envisage that the effect of the thermal fluctu- 
ations in the present case due to heating is not as big as 
when rj = to overcome the correlation in extent that it 
reduces the degree of squeezing visibily. 

It goes beyond contention that the degree of squeez- 
ing in such a system is substantially limited by the ther- 
mal fluctuations specially at earlier stages and relaxed 
somewhat as time progresses since the strength of the 
correlation gets better as more atoms contribute towards 
establishing the correlation. Predominantcly, it is the 
competition between the effect of the thermal fluctua- 
tions and the nonclassical correlation of the two modes 
that lead to the characteristic evolution of the two-mode 
squeezing and other nonclassical features. 
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In order to understand the existing situation indepth, 
one needs to revisit the master equation (ffU|). From the 
general form of the master equation, it is not difficult to 
infer that the second term stands for the gain of mode a, 
the third term for the lose of mode b and the first term 
for the lose of mode a. It is also noteworthy that the first 
term and the part involving k in the third term comes 
from the contribution of the coupling of the cavity with 
external vacuum modes (reservoir) and basically related 
to the thermal fluctuations. On account of this, it may 
not be difficult to observe that until the lasing system 
able to generate strong light, the effect of the thermal 
fluctuations dominates and later, undoubtedly, the effect 
of the lasing system begins to overtake it. The result ob- 
tained here, the two-mode squeezing would be maximum 
at steady state, corroborates with this line of reasoning. 

It is vividly presented in Figs. [2] and [3] that the de- 
gree of two-mode squeezing is significantly impaired by 
chaotic light, that is, the more intense the seeding, the 
more the squeezing is damaged. It turns out that there is 
no essential difference between when the initial thermal 
light is in mode a and mode b. Referring to Eqs. 
(Unj) and ([50)1 bespeaks that at t = 0, the thermal seed 
contributes to the variances of the quadrature operators 
(Acj_) a term equal to n a + h^. In other words, at t = 0, 
Ac 2 . = I + n a + ni,. That is why, in the vicinity of t = 0, 
the generated radiation would be in squeezed state only 
when the strength of the seed is significantly weaker. 

In the later stages of the process, the correlation of the 
emitted radiation would be strong enough to overcome 
the effect of the seed. If one critically observes Figs. [H [3] 
and[4l one may see that the degree of two- mode squeezing 
is almost the same for all cases after t = 2. This suggests 
that the effect of the seed is relatively short lived. 

Furthermore, making use of Eqs. (gS]), (gHJ), ([50]) and 
(|52|) . one readily finds at steady state 

Ac| =1 + ^(1-7?) 

|~4k(k + An) + A[A + (2k + Ar])](l ± y/l-r, 2 ) 
4k(2k + At))(k + Aif) 

± Ay/ 1 - rf 

x 4k(k + Arj) +A 2 {r ] 2 ~lTVl~ V 2 ) ] ^ 
4k(2k + Arj)(n + A?]) 

It is not difficult to realize from Figs. [21 El HI and Eq. 
(f53|) that the degree of squeezing at large time scale is 
independent of the mean photon number of the chaotic 
light. This is mainly because as time progresses there are 
more atoms that pass through the cavity in a way absorb 
the available radiation. Since the number of chaotic light 
in the cavity cannot increase, the absorption mechanism 
dwindles it over time. 

At this juncture, it would be appropriate noting that 
the atoms are continuously injected with constant rate 
and left the cavity after decaying to any other energy 
level that does not involve in the lasing process. In con- 
nection to this, what most probably happens is an atom 



after absorbing a chaotic light emits a light which is not 
necessarily has the same nature as it absorbs. Since this 
process successively repeated, the chaotic light in the cav- 
ity is converted to a correlated light or lost all the way 
which leads to diminishing effect of the seed on the prop- 
erties of the cavity radiation. 

B. Mean photon number 

The mean number of photon pairs describing the two- 
mode cavity radiation can be defined as 

N = {<*(t)c(t)), (54) 

where c(t) is the annihilation operator (|3Tj) . It is read- 
ily noticeable that the operators in Eq. (|54p are in the 
normal order. Hence it is possible to express Eq. (|54[) in 
terms of c-number variables associated with the normal 
ordering as 

N=±[{a*(t)a(t)) + {p*(t)m)], (55) 

where 2N represents the mean photon number 

In the following, the mean photon number is plotted 
against time for rj sa 0, k — 0.5 and different values of A, 
n a and fib- Here the value of rj is chosen on the basis of 
earlier studies at steady state in which the mean photon 
number is found to be maximum for rj = 0. The value 
of A on the other hand is fixed to be small so that the 
dynamics of the intensity of the radiation be evident from 
the plots. 
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FIG. 5: Plots of the mean photon number of the cavity ra- 
diation (2JV) for k = 0.5, n a — fit — 0, r\ ~ 0, and different 
values of A. 

It is clearly shown in Fig. [5] that the mean photon 
number increases with time and linear gain coefficient 
when there is no seeding. This is readily understandable 
since traversing more atoms across the cavity most likely 
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FIG. 6: Plots of the mean photon number of the cavity radi- 
ation (2JV) for k — 0.5, fib = 0, f] w 0, A = 1, and different 
values of n a . 



FIG. 8: Plots of the mean photon number of the cavity radi- 
ation (27V) for k — 0.5, r\ m 0, A = 1, and different values of 
n a = fib- 
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FIG. 7: Plots of the mean photon number of the cavity radi- 
ation (2N) for k = 0.5, n a = 0, r\ » 0, A = 1, and different 
values of n;,. 



leads to generation of more photons. However, the results 
indicated in Figs. [H [7J and [5] do not support this idea 
when there is a seed in general terms. 

To begin with, when fib = and h a ^ 0, the situation 
indicated in Fig. |S1 the mean photon number increases 
with time at the beginning and then decreases later. This 
is mainly related to the fact that at initial time the atoms 
emit until the radiation in the cavity is sufficient enough 
to frustrate the transition from energy level \a) to \b). 
Later due to the availability of large number of photons 
than the system can generate, the atoms begin to pickup 
these photons leisurely and underwent emissions to any 
other level that does not involve in the lasing mechanism. 
In the process, as earlier results display, the number of 
chaotic light starts to dwindle which compels the decline 



of the mean photon number. Critical scrutiny reveals 
that the mean photon number does not slump to the 
saturation value by more than the mean photon number 
initially in the cavity and the mean photon number can- 
not be less than its value at t = for parameters under 
consideration. 

Quite contrary to this, when n a = and n& ^ 0, the 
situation dipcited in Fig. [71 the mean photon number ini- 
tially decreases and later starts to increase with time. In 
this case, it is worth noting that nearly 50% of the atoms 
are prepared to be in the lower energy level and most 
probably hungry for photons. As a result, the photons 
due to the seed which are available in the cavity read- 
ily absorbed by this atoms and consequently the number 
of the photon drastically falls. As critical scrutiny might 
have revealed, the depletion in mean photon number can- 
not go lower than the initial mean photon number for all 
conceivable cases. 

After quite sometime, the available photons (fib) be- 
gins to dwindle whereby the usual emission process gets 
a chance to regain control in which the mean photon 
number readily increases till it becomes nearly uniform 
at longer time scale. It is not difficult to observe from 
Fig. [5] that specially when n a = fib = 1.5, both these ef- 
fects are readily observed. That means the mean photon 
number decreases with time, then increases and later de- 
creases with time. This indicates that the absorption of 
mode b is efficient in earlier times whereas that of mode 
a is at later times. Almost similar behavior of mean pho- 
ton number has been predicted for different parameters 
when the initial radiation is assumed to be a coherent 
light [|. 
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V. QUANTIFICATION OF ENTANGLEMENT 

Based on different physicai contexts and mathemati- 
cal considerations, several sufficient inseparabilit y c rite- 
ria for composite state have been proposed |20M3G | for 
continuous variables. One of these alternatives is the cri- 
terion associated with the logarithmic negativity [311433T ] 
that can be defined as 



Furthermore, making use of Eqs. ([58]). (|59j) . and the 
accompanying definitions of m, n and c, one can readily 
see that 

detA = 4(a*(t)a(t))[(a*(t)a{t)) + 1] + 1, (61) 



detB = 4(P*(t)/3(t))l(P*(t)Kt)) + 1] + 1, 



(62) 



E N = max[0, - log 2 V s ], 



(56) 



where V s is the smallest eigenvalue [2CJ, |31[ . 

Upon solving the eigenvalue Eqs. for symplectic spec- 
trum of the covariance matrix of the partially transposed 
density operator, the smallest eigenvalue is found to have 
the form 



V s = 



a - \Ja 2 - idetfl 



-1 1/2 



(57) 



with a = detA + detB — 2detCAB, in which A and B 
are the covariance matrices that represent each modes 
separately while Cab is the intermodal correlations and 
ft is the covariance matrix to be defined later [gTj, H3] • 

In line with this consideration, the 2X2 block form of 
the covariance matrix can be expressed as 



= 



A C AB 
B 



^AB 



(58) 



where Sly = ^(X^Xj + XjXj) — (Xi)(Xj). It is essential 
to note that the involved quadrature operators can be 
defined as X\ = a + a\ X2 = i \a> — a] , X3 — b + b' and 
X4 = i\b' — b]. It is also worth noting that Xj's can be 
related to the quadrature operators by c+ = -^(Xi+Xs) 
and c_ = -^(X 2 +X 4 ). 

With this introduction, the extended covariance ma- 
trix turns out to have the form 



n 



m c 

m — c 

c n 

-c n 



(59) 



where m = 2(a t (t)a(t)) + 1, n = 2(S t (£)6(£)) + 1 and 
c = (d(t)b(t)) + (&(t)b'(t)) which can also be expressed 
in terms of c-number variables associated with the normal 
ordering as m = 2(a*(t)a(t)) +l,n= 2{(3* (t)(3{t)) + 1 
and c = (a(t)(3{t)) + (a* (t) /3* (t)) . 

The two-mode Gaussian composite state is entangled 
when the logarithmic negativity is positive, that is, if 
En = —Log?y s which implies that Log^Vg should be 
negative. This, on the other hand, holds true provided 
[20t that 



V s < 1. 



(60) 



detCAB 



-4(a(t)(3(t)) 2 , 



(63) 



detfl = 16 



(a*(t)a(t)> + </3* (t)/3(*)) 



(a*(t)a(t))W* (t)P(t)) ~ Ht)m? 



(64) 



In order to study the evolution of the entanglement us- 
ing criterion (|60|) . the smallest eigenvalue (V s ) is plotted 
against time applying Eqs. (g5J|, (|5U|) . (|6Tj) . 

dill), and ([54]) . Taking the value of 77 to be very 
close to zero and similar values for the rest of the pa- 
rameters as for two-mode squeezing, the evolution of the 
entanglement would be studied. 




>" 0.7 



FIG. 9: Plots of the minimum of the logarithmic negativity 
(V s ) for k = 0.5, 77 ~ 0, h a — fib = 0, and different values of 
A. 

It is clearly indicated in Fig. [3] that the degree of en- 
tanglement increases with the linear gain coefficient for 
values under consideration. It so happens that the de- 
gree at which the criterion (|60]) is violated increases with 
time at early stages of the operation, but it decreases 
in the later stages. Similar behavior of the nonclassical 
features including two-mode squeezing fl8l ] and entan- 
glement measure following from Duan et al. [j| Hl[ when 
there is an external driving radiation and phase fluctua- 
tions was reported. 

Most recently, a thorough explanation in this regard 
is provided in relation to the thermal fluctuations due 
to vibrations of the atoms on the walls of the cavity. 
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FIG. 10: Plots of the minimum of the logarithmic negativity 
(V s ) for k = 0.5, 77 ss 0, fib = 0, A — 10, and different values 
of n a . 



FIG. 12: Plots of the minimum of the logarithmic negativity 
(V s ) for n = 0.5, r\ « 0, A — 10, and different values of 
n a = fib- 
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FIG. 11: Plots of the minimum of the logarithmic negativity 
(V s ) for k = 0.5, T) » 0, n a = 0, A = 10, and different values 
of fib. 



Complementary explanation is also provided in Section 
IV. Though it is not clearly observed here, earlier stud- 
ies show that the minimum of the logarithmic negativ- 
ity can go as small as 0.25 for A = 10 at steady state 
PJj. Moreover, comparison of Figs. [5] and M reveals 
that V s attains minimum value (better degree of entan- 
glement) quite earlier than when the intensity becomes 
strong. Therefore, apart from what is expected at steady 
state, as it stands, the practical utility of this system for 
generation of a strong entangled light does not seem sat- 
isfactory. 

Furthermore, it is not difficult to see from Figs. [10] 
and [Tl] that at the vicinity of t = 0, the entanglement 
of the cavity radiation when n a = 0; Sj, ^ and 
7% a 7^ 0; fib — exhibits different nature. By large, 



this result can be associated with the nature of the mean 
photon number in a similar case. Even though the pres- 
ence of thermal seed in mode a somewhat diminishes the 
degree of entanglement, its effect would be more promi- 
nent when the thermal light is seeded in mode b. Critical 
comparison between Figs. and [TU1 for A = 10, indi- 
cates that the effect of the thermal light when it is in 
mode a is quite nominal over time, but as clearly shown 
in Fig. [T2]the effect would be large when n a — fib is large 
specially at the beginning of the lasing process. In conso- 
nant with earlier observations, the effect of the thermal 
seeding would be insignificant in the later time of the 
operation to change the achievable entanglement. 

In order to see the relation between two-mode squeez- 
ing and entanglement, V s and Ac 1 ! are plotted against 
time for selected parameters. It is evident from Fig. [T3] 
that the nature of the evolution of the degree of entan- 
glement and two-mode squeezing is basically similar, al- 
though observing entangled light when there is no squeez- 
ing at early stages of the process is possible. The resem- 
blance as in Fig. [13] remains appearing at larger time 
scale as predicted earlier at steady state [13| . It is worth 
noting that such a tendency of exhibiting similar behav- 
ior would break near rj = 0, since the squeezing disap- 
pears at rj = while the degree of entanglement is nearly 
close to 50% which is undoubtedly a significant difference. 
It is due to this disparity that the values of rj are taken to 
be different while investigating two-mode squeezing and 
entanglement. 



VI. CONCLUSION 

In this contribution, detailed analysis of the time evo- 
lution of the cavity radiation of the two-photon correlated 
emission laser when the cavity is initially maintained at 



11 



2.5 




2 4 6 8 10 

time 

FIG. 13: Plots of the minimum of the logarithmic negativity 
(Vs) and the minus quadrature variance (Ac?.) for k = 0.5, 
r\ — 0.2, A = 10, and n a = n(,=0.5. 

a particular temprature is presented. The effect of the 
seed thermal light is included via the parameters that 
describe the initial conditions in the solution of the cor- 
responding stochastic differential equations. With the 
understanding that the effect of the thermal light enters 
as initial conditions only if the source is entirely shut off 
prior to the injection of the atoms, the characteristic cor- 
relations of the chaotic light are employed in the process 
of determining various quantities of interest. 

It turns out that the thermal light significantly dam- 
ages the nonclassical features like two-mode squeezing 
and entanglement at earlier stages of the lasing pro- 
cess. Although the degree of entanglement and two-mode 
squeezing exhibit more or less a similar dependence on 
the intensity of the thermal light and time, their nature 
comes out to be different when the thermal light is seeded 
only in one of the modes specially at the vicinity of t = 0. 
Fundamentally, the degree of two-mode squeezing is al- 
most independent of which mode is initially seeded, but 



the degree of entanglement reduces considerably when 
the light with the same strength is seeded in mode b. 

Quite remarkably, a somewhat similar characteristic is 
perceived for the mean photon number. When viewed 
over a longer time span, the mean photon number has 
a peak when mode a is seeded, but a dip when mode b 
is seeded. The time for which one observes the peak is 
relatively much greater than the corresponding time for 
the dip. It is found that this disparity in the tendency of 
the mean photon number increases with the intensity of 
the seed. 

Moreover, it is unambiguousely shown that at longer 
time scale, the two-mode squeezing, entanglement and 
intensity of the radiation are independent of the strength 
of the initial thermal light. This by large is related to the 
presumption that the atoms can absorb any light in the 
cavity as long as it is resonant. In view of this, it would be 
an obvious matter to realize that the effect of the thermal 
light is sucked from the cavity in a level that it would not 
affect the nonclassical features of the radiation. 

As long as the radiation in the cavity is resonant with 
the radiative levels of the atoms, it is not difficult to infer 
that this outcome remains valid irrespective of the type 
of the seed. This entails that while designing the exper- 
imental setup, one should not be compelled by the rem- 
nant radiation in the cavity, since its effect would be sig- 
nificantly nominal at longer time scales. It is hoped that 
such an outcome can provide another dimension to the 
practical utilization of the various schemes of the three- 
level laser as a source of bright entangled light. 
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